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ABSTRACT 



We calculate all point symmetries of the Fokker - Planck equation in two-dimen- 
sional Euclidean space. General expression of symmetry group action on arbitrary solu- 
tion of Fokker - Planck equation is presented. 



1. The symmetries of the Fokker - Planck equation in two dimensions 

The object of our considerations is a special case of Fokker - Planck equation, which describes evolu- 
tion of 2D continuum of non-interacting particles imbedded in a dense medium without outer forces. The 
interaction between particles and medium causes combined diffusion in physical space and velocities space. 
The only force, which acts on particles, is damping force proportional to velocity. 

The 3D variant of this equation was investigated in our work [1]. In this work fundamental solution 
of 3D equation was obtained by means of Fourier transform. 

The ID variant of this equation was investigated in our work [2]. All point symmetries of the Fokker 
- Planck equation in one-dimensional Euclidean space were calculated. 

In present work we continue this investigation for more complex 2D equation. 

The Fokker - Planck equation in two dimensions is 



where 

n — n{t, X, y, u, v) - density; 

t - time variable; 

x,y - space coordinates; 

M, V - velocity; 

a - coefficient of damping; 

k - coefficient of diffiision. 

The list of symmetries of the Fokker - Planck equation in one dimension follows. The calculations of 
symmetries are rather awkward. They are carried out to APPENDIX 1 . 

Instead of classic - (j)" notation we use another ("^") notation. This notation was presented in our 
work [2]. 

Addition of arbitrary solution 




(1) 



dn ' 



(2) 
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where A is arbitrary solution of the (1) equation. 
Scaling of density 



\2 = n 



dn 



The reason of symmetries (2-3) existence is linearity of PDE (1). 
Time shift 



V3 = 



dt ' 



Space translations 



Space rotation 



^ du " 3v dx ^ dy 

Transformations (5) and (6) build two-dimensional EucUdean movements group. 

Extended GaUlean transformations, which besides time and space coordinates affect the density 

d d an d d d an 3 

V7 = — + f (ax + u)^; Vg= — + f^; (ax + v) 

on 



du dx 2k 



dx 2k 



dn 



(3) 



(4) 



(5) 



(6) 



(7) 



Negative exponent transformations - they affect time and space coordinates, contain time-dependent com- 
mon multipUer (negative exponent). They do not affect density. 



V9 = e 



du dx / ^ { ^ dv dy J 



(8) 



Positive exponent transformations - they affect time, space and density, contain time-dependent common 
multiplier (positive exponent). 

d^ J d d _a^ d 

dv dy k dn 



Vn - e"" \ a + - — nu ^ 1; V12 - e"'] a i^— + — mv 

du dx k dn 



(9) 



One-parameter groups, generated by vector fields Vj - V12, are enumerated in the following list. The 
list contains images of the point (n, t, x, y, u, v) by transformation expieVj) 



Gl: 

Gy. 
G4: 



{n + eA,t,x, y,u, v); 
(e^n, t, X, y, u, v); 
{n,t + e, X, y,u, v); 
{n,t, X + £,y,u, v); 
(n, t,x,y + e, u, v); 

(n, t, cos(£)x + sin(£)>', - sin(£)x 4- co^ie^y, cos(£)m -I- sin(£)v, - sin(f )m -I- cos(£)v); (10) 
1 



(exp 



^ e(ax + m) + ^ £^(a/ + 1) 



n,t,x + et, y,u + e, v); 
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G9: 
GlO: 

Gil: 
Gi2: 



(exp 



n, t,x,y + et, u,v + e); 



(n, t,x + ee "' , y,u - ea e "', v); 
(n, t,x,y + ee'"', u,v — ea e~"'y, 
1 



(n exp 



(n exp 



a 



— ei eu+^e^a e"' 



-e"'\ ev + -e^a e"' 
k { 2 



,t,x + ee"', y,u + ea e"', v). 



,t,x,y + ee" ,u,v + ea e" ). 



For relatively nontrivial integration of Gj ,Gg , Gn and G12 we refer to [2]. 

The fact, that G, are symmetries of PDE (1) means, that if f{t, x, y, u, v) is arbitrary solution of (1), 
the functions 



,,(1). 



,,(2). 



.,(8). 



.,(12). 



fit, X, y, u, v) + eA{t, x, y, u, v); 
e^f{t,x, y,u, v); 
f{t — e, X, y, u, v); 
f{t, X — e,y, u, v); 
fit, x,y-e, u, v); 

f{t, co%{e)x - %m{e)y, %\a{e)x + cos(£)>', cos(£)m - sin(£)v, sin(£)M + cos(£)v); 
1 



exp 



exp 



e{ax + u)--e^{at+l) 



2^ ^e{ay + v)--e\at+\) 



f{t, X - et,y,u - e, v); 



f{t,x,y-et,u,v-ey. 



(11) 



fit, X - ee'"', y,u + ea e'"', v); 
/(/, x,y- ee~"', u,v + ea e'"'); 
1 



exp 



exp 



e"' \ £u e^a e"' 

k I 2 



~e'"\ev- -e''a e'' 
k [ 2 



fit, X - ee" ,y,u- eae" , v). 



/(/, x,y- ee"' , u,v- eae"'). 



where e - arbitrary real number , also are solutions of (1). Here A is another arbitrary solution of (1). 

We systematically replaced "old coordinates" by their expressions through "new coordinates". Note, 
that due to these replacements terms with e^ in m^^' , u^^^ , m^'^ and m^'^' change their signs. 
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We have trivial solution n - e^"' at our disposal. If we act on this solution by transformations (10), 
we obtain 4 new solutions: 



n - exp 



n - exp 



2at - — I e{ca + m) - - E'^{at + 1) 



2ar - ^ \e{ay + v) - ^ e^{at + 1) 



n — exp 



n — exp 



(12) 



(13) 



(14) 



(15) 



General expression is 

U = e"^ exp 



exp 



of _ _ I 2 

— I £:8(a3; + v)- -£8(a/+ 1) 



X (16) 



xexp 



f 1 
-'^e'"ieniu-e7)--siiae"' 



exp 



xf{t - e^, X - £4 - Eft — egc - Sn^"', y-£5- egt- - Sne"', m - £7 + e^a e - Sua e"', v-eg + ei^a e - e^a e"') + 

+£iA{t, X, y, u, v); 

where 

X = cos(£6)x - sin(fg)>'; (17) 

y = sin(£6)x + cos(£6)y; (18) 

u - cos{e^)u - sin(£6)v; (19) 

V - sin(£6)M + cos(£6)v. (20) 
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DISCUSSION 

Looking at the list of all point symmetries of the Fokker - Planck equation in two-dimensional Euclidean 
space, we see, that there is no simple way to get, for example, fundamental solution of PDE, using these 
symmetries. We have not at our disposal such an instrument, as scaling of independent variables t, x, y, u, v. 
The result (12-15) of action of symmetry group on trivial solution is not very interesting from physical 
point of view. 

Indirect way of use of Galilean transformations (7) was demonstrated in [1]. The transformation was 
used for generalisation of solution, which was obtained in the form of exponent of quadratic form of space 
coordinates and velocities with time dependent coefficients. 

There is need of further investigations of Fokker - Planck equation and its set of symmetries, which 
may lead to another physically interesting results. We can follow the scheme of [7J ; to consider invariant 
solutions for some one-parameter group, thus reduce the independent variables number. To find for 
obtained in such a way equation all point symmetries - and so long. 

In the work [7] this scheme was represented for equations of elasticity and plasticity. 
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APPENDIX 1 



The infinitesimal invariance criteria for PDE (1) is 

dn dn dn dn dn 

S{^) + Su^ + uS{^) + Sv ^ + vd(—)- (Al-1) 
at ax ax ay dy 

dn dn dn dn d^n d^n 

—aSu ^ auS{^) — aSv avS{^) — laSn - k{S{^r-^) + S{^r-^)) - 0. 

du du dv dv du^ dv^ 



According to [2] (APPENDIX 1, eq. (Al-8) and (Al-18)), we have for variations of derivatives fol- 
lowing expression: 

dn d d dn dn d d dn 

<J — = — {Sn) + — {Sn) ^ - ^ {5x) + — (Sx) — ) - (Al-2) 
ax ax an ax ax ax an ax 

dn d d dn dn d d dn 

- ^ {5y) + — {5y) ^) - ^ {5u) + — {5u) — ) - 
dy dx dn dx du dx dn dx 

dn d d dn dn d d dn 

- ^ iSv) + — (^v) ^) - ^ (St) + ^ (St) — ); 
dv dx dn dx dt dx dn dx 

dn d d dn dn d d dn 

S^^^ (Sn) + ^ (Sn) ^-^{^ (Sx) + ^ (Sx) ^) - (Al-3) 
dy dy dn dy dx dy dn dy 

dn d d dn dn d d dn 

dy dy dn dy du dy dn dy 

dn d d dn dn d d dn 

dv dy dn dy dt dy dn dy 

dn d d dn dn d d dn 

S^^^ (Sn) + ^ (Sn) ^ - ^ (Sx) + ^ (Sx) ^) - (Al-4) 
du du dn du dx du dn du 

dn d d dn dn d d dn 

dy du dn du du du dn du 

dn d d dn dn d d dn 

dv du dn du dt du dn du 

dn d d dn dn d d dn 

S^ = ^ (Sn) + ^ (Sn) 5" - 3" (5- (Sx) + ^ (Sx) 3-) - (Al-5) 
dv dv dn dv dx dv dn dv 

dn d d dn dn d d dn 

ay av an av au av an av 

dn d d dn dn d d dn 

- — (- {Sv) + — {dv) ^) - ^ {St) + — {St) —); 
av av an av at av an av 

dn d d dn dn d d dn 

^ ^ = ^ {Sn) + — {Sn) ^ - ^ {Sx) + ^ {Sx) —) - (Al-6) 
at at an at ax at an at 

dn d d dn dn d d dn 

dy dt dn dt du dt dn dt 
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dn d d dn dn d d dn 

dv dt dn dt dt dt dn dt 



,3^ 



dn du^ dudx du 



d^n d d dn 

{—(dx)+ — {Sx)—) 

an au auay au 



d^n d d dn 

{—{5y)+ — {5y)—)- (Al-7) 
an au 



d^n d d dn 

du?- du dn du dudv ^du 



d^n d d dn 

dn du dtdu ^3m 



d^n d d dn 

i—m+^m^)- 

an au 



d^n dn d d^n dn d d^n dn d d^n dn d d^n dn d 

du^ dx dn du^ dy dn du^ du dn du^ dv dn du^ dt dn 

d^ d^ dn dn d^ d^ dn 

+ ^ (Sn) + 3-3- (Sn) ^ + 3- (5-3- (Sn) + ^ (Sn) ^) - 
au-^ anau au au anau an^ au 



dn d^ 
dx du^ dndu 



dn 

(Sx) 3-)- 
au 



dn 32 32 3^ 

dy du^ dndu du 



^(^^^Su) + — 
du du^ dndu 



dn 
au 



dn d^ d^ dn 

dv du^ dndu du 



dn ^ d^ ^^^^ ^ d^ 
dt du^ dndu 



dn 
du 



dn dn d^ d^ dn 

JT JT (fit + ^ ir^ - 

dx du dndu dn-^ du 



dn dn ^ d^ (S ) (S ) '^") 

dy du dndu dri^ du du du dndu 



dn dn d^ d^ dn 

{^^{5u)+ — {du) — )- 
an-^ au 



dn dn d^ d^ dn 

dv du dndu dn^ du dt du ^dndu 



dn dn d^ d^ dn 

{.^{5t)+^{5t)^)- 
dn^ du 



d^n d d dn d^n d d dn 

■ J— (j- (Sx) + — (Sx) — ) - 3— (— (dy) + ^ (Sy) j-) - 
auax au an au auay au an au 



d^n d d dn 

du^ du dn du dudv ^du 



d^n d d dn 

dn du dtdu ^du 



d^n d d dn 

dn du 



3^ 
dv^ 



d^n 



dn ^ dv^ dvdx 



d^n ^ d 



d dn 
{Sx)+ — {Sx) —)- 
dn dv 



d^n d d dn 

^i^iSy)+^iSy)^)- (Al-8) 
dvdy dv dn dv 



d^n d d dn d^n d d dn 

dudv dv dn dv dv^ dv dn dv dtdv ^dv 



d^n d d dn 

dn dv 



d^n dn d 
dv^ dx dn 



(Sx)- 



d^n dn d 
dv^ dy dn 



(Sy) - 



d^n dn d 
dv^ du dn 



(Su) - 



d^n dn d 
dv-^ dv dn 



(dv)- 



d^n dn d 
dv^ dt dn 



(5t) + 



d^ d^ dn dn d^ d^ dn 

+ J— (Sn) + J— (Sn) j^ + jr- ( viT + ^ JT^ ' 
dv^ dndv dv dv dndv dn^ dv 



dn 32 32 dn dn 3^ 3^ 

dx dv-^ dndv dv dy dv^ dndv 



dn 



_dn d^ 32 
du dv^ 3n3v 



dn 
(Su)^)- 



dn 32 32 dn 

dv dv^ dndv dv 
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dn dn dn dn dn 

dt dv^ dndv dv dx dv dndv dn^ dv 

dn dn 9^ 3^ dn dn dn d^ d^ dn 

ay av anav an^ av au dv dndv dn^ dv 

dn dn ^ d^ (S )+ (S ) ^" ^" ( (St)+ {St)^^) 

dv dv dndv dri^ dv dt dv dndv dn^ dv 

d^n d d dn d^n d d dn 

dvdx dv dn dv dvdy dv dn dv 

d^n d d dn d^n d d dn d^n d d dn 

dudv dv dn dv dv^ dv dn dv dtdv dv dn dv 

dn 

We eliminate ^ in (Al-1) using original equation 
dt 

dn ( dn dn dn dn d'^n d^n ^ 

Collecting similar tenns, we obtain following equations: 
dn dn 
dx du 

d^ d' 

-2ku^^^{5t) + 2k^^^{dx)^0; (Al-10) 
dndu dndu 

dn dn 
dx du^ 

-ku^iSt) + k^iSx)^0; (Al-11) 
dn^ dn^ 

dn dn 
dx dv 

-2ku (St) + 2k (dx) = 0; (Al-12) 
dndv dndv 



dn dn 
dx dv^ 



dn 
dx 



^2 ^2 

-ku^(St) + k^(Sx) = 0; (Al-13) 
dn-^ dn^ 



d d 3,3 3 

-auv — (St) + 2aun — (St) + au — (Sx) - au — (St) + av — (Sx) - (Al-14) 
dv dn du du dv 

3 32 32 32 32 

-2an — (Sx) - ku (St) - ku (St) + k (Sx) + k ^r-^ (Sx) + 
dn du-^ dv^ du'^ dv^ 

d d 3,3 3 3 

+MV ^ (^0 - M ^ (Sx) + M ^ (^0 + ^ (^0 - V — (Sx) + Su- — (Sx) ^ 0; 
dy dx dt dx dy dt 
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dn d^n 



du dudx 



dn d^n 



du dudy 



dn d^n 



du dtdu 



d 

2k ^ (Sv) = 0; 
dn 



d 

2k — (Sx) = 0; 
dn 



2k^iSy)^0; 
dn 



d 

2k — {St) = 0; 
dn 



dn 
du 

3^ 3^ 3^ d 

aku ^r-^ {St) + aku ^r-^ {St) + Aakn ^ ^ {St) — auv ^{St) + 
du^ dv^ dndu dy 

d 3,3 3 3 

+au -zr- {Su) — au^ {St) — au^ {St) + av ^ {Su) — 2an ^ {Su) — aSu + a 
du dt dx dv dn 

,3 , , 3 3^ 3^ 3^ 3 

-2a^un — {St) + a^u^ — {St) - 2k {Sn) + k — {Su) + k — {Su) -u — {Su) - \ 
dn du dndu du^ dv^ dx 

dn dn 
du^ dv 

akv — ^ {St) + k^ {Sv) = 0; 

dn d^n 
du^ du?- 

2 3' 

dn d^n 
3m2 3v2 

2 3' 

^^a;^2(^0^0; 

dn 
d^ 

^2 ^2 ^2 g2 

2aku {St) + 2akn ^ {St) -k^ {Sn) + 2k {Su) = 
dndu dn^ dn^ dndu 

dn 
du^ 

32 32 
aku zr-^ {St) + k zr-^ {Su) - 0; 
dn-^ dn^ 



-11- 



dn d^n 



dn d^n 



dn d^n 



9v dudv 



dn d^n 



dv dvdx 



dn d^n 
dv dvdy 



dn d^n 



dv dtdv 



dn 
dv 



(St) + (St) 0; (Al-34) 

onav anov 



d d^ 

2k^(Sv) + 2k^ ,-^(St) = Q; (Al-35) 



-av 



dn dndv 



d 

2k ^ (Su) ^ 0; (Al-36) 
dn 



d 

2k — (Sx) = 0; (Al-37) 
dn 



d 

2k ^ (Sy) = 0; (Al-38) 
dn 



d 

2k — (St) ^ 0; (Al-39) 
dn 



d^ d^ d^ 

akv (St) + akv (St) + 4akn :-— (St) - (Al-40) 
dM^ dv^ dndv 

3 3 3 3 

—auv zr- (St) + au ^ (Sv) + av zr- (Sv) — av^r- (St) — 
dx du dv dt 

2 ^ (St) - 2an ^ (Sv) - aSv + a^uv ^ (St) - 2a\n ^ (St) + 
dy dn du dn 



+a\'^ ^ (St) + (Sv) - 2k —— (Sn) + k^ (Sv) - m — (Sv) - v — (Sv) - ^ (Sv) = 0; 
dv dw^ dndv dv^ dx dy dt 



dn d^n 
3v2 3m2 



dn d^n 
3v2 dv^ 



g2 

^ (St) = 0; (Al-41) 



g2 

k^ (^0 0; (Al-42) 



3n2 
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32 32 32 

2akv (St) + lakn ^ {St) -k — (Sn) + 2k z-z- (Sv) 
dndv an^ dn^ dndv 



3^ 3^ 
akv ^ (St) + k^ (Sv) = 0; 



3 3 3 3 

aku — (St) + akv — (St) - lakn — (St) -ku — (St) 
(III dv an ax 

3 3 3 3*^ 3^ 

-kv ^ (St) + 2kz- (Su) -k^ (St) + k^^ (St) + k^ ^ (St) ■■ 
ay au at au^ av^ 



3 3 3 3 

aku ^ (St) + akv -zr- (St) — lakn zr- (St) -ku — (St) — 
au av on dx 

3 3 3 ,32 ,32 

kv — (St) + 2k — (Sv) -k^ (St) + /t^ — (St) + k^ — (St) 
dy dv at du^ dv^ 



3 3 
2k ^ (Su) + 2k^ (Sv) = 0; 
av au 



2k ^ (Sx) = 0; 
au 



2k^(Sy) = 0- 
du 



3 

2k ^ (Sx) = 0; 
dv 



2k^(Sy)^0; 



-13- 



dtdv 



d 

2k ^ (St) = 0; (Al-52) 



d 

2k — (St) = 0; (Al-53) 
av 



2akn ^r-— (St) + 2akn ^r— (St) - 2aun -z— (St) — au^ (Sn) — (Al-54) 
3m^ av^ ax au 

3 3 3 3 

—2avn :r— (St) — av — (Sn) + 2an zr- (Sn) — 2an z— (St) — 
dy dv dn dt 

-2aSn + 2a^un ^ (St) + 2a^vn ^ (St) - Aa^n^ ^ (St) -k^ (Sn) - 
du dv dn du^ 

3^ 3 3 3 

-k z—r (Sn) + M (Sn) + V (Sn) + (Sn) = 0. 
av^ dx dy dt 

From (Al-37 - Al-39), (Al-48 - Al-51) we see, that Sx = Sx(x,y,t); Sy = Sy(x,y,t); 
St - St(x, y, t). Using these expressions, we simphfy the rest of equations (Al-10 - Al-54). 

3 

I + - ^ (Sx) = 0; (Al-55) 
dt 



3 




3 




3 2 


3 3 


dy 


(St)- 


dx 


(Sx) + u 


^^(St)^u^ 


^(St)-v^< 
dx dy 


3 




3 




3 


3 , 3 


3a: 


(St)- 


dx 


(Sy)-v 


^(Sy) + v 
dy 


3^^"^^^^^ 3^ 



3 

dt' 



(Al-56) 



-auv ^ (St) + au z— (Su)- au zr- (St)- au zr-(St)+ (Al-57) 
dy du dt dx 

3 3^ 3^ 

+av zr- (Su) - aSu - 2k .. .. (Sn) + k zr-ir (Su) + 
dv dndu du^ 

, 3^ 3 3 3 

+k zr-^ (Su) — Uzr—(Su) — Vzr- (Su) — zr- (Su) = 0; 

dv^ dx dy dt 

32 

-k ^ (Sn) = 0; (Al-58) 
dn^ 

3 3 3 3 

-auv — (St) + au — (Sv) + av — (Sv) - av — (St) - (Al-59) 
dx du dv dt 



-aV^ rr- (St) -aSv + kr—r (Sv) " 2k (Sn) + 

dy du-^ dndv 

32 3 3 3 

+k ^ (^v) - « ^ (^v) - V ^ (Sv) - ^ (Sv) = 0; 

32 

-k ^ (Sn) = 0; (Al-60) 
dn^ 
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o 




d d 


a 
o 


dx 




-kv^ (dt) + 2k^ (du) - 
oy ou 


at 


d 

dx 


(St) 


d d 

-kv^ (St) + 2k^ (Sv) - 
dy dv 


d 
dt 






d d 








2k — (Su) + 2k — (Sv) = i 
dv du 


3; 



1 = 0; (Al-61) 
I = 0; (Al-62) 
(Al-63) 

a a a a 

-2aun -zr- (St) -au zr- (Sn) - 2avn ^ (St) -av ^ (Sn) + (Al-64) 
dx du dy dv 

a a a^ 

+2an ^ (Sn) - 2an ^ (St) - 2aSn -k^ (Sn) - 
dn dt du^ 

a^ a a a 

-k r- (Sn) + U,— (Sn) + v — (Sn) + — (Sn) = 0. 
dv-^ dx dy dt 

From (Al-58) and (Al-60) we conclude, that 

Sn = A + nB; (Al-65) 

where A = A(x, y, u, v, t), B - B(x, y, u, v, t). Using this expression, we simpUfy the rest of equations 
(Al-55 - Al-64). 



a 




a 




a 2 


a a 


d 






a^ 


(Sx) + u 


^^(St)^u^ 


^ (St) - V ^ (Sx) + Su- 
dx dy 


dt 


a 




a 




a 




a 


a^ 


(St)- 


dx 


(Sy)-v 


^(Sy) + v 
dy 


^^(St) + v^^(St) + Sv- 


dt 



-auv -zr- (St) + auzr- (Su) -auzr- (St) -au zr— (St) + (Al-68) 
dy du dt dx 

a a^ a^ as a a a 

+av zr- (Su) — aSu + k zr-^ (Su) + k ^r-— {Su) — 2k Uz— {Su) — v^r- (Su) — ^ (Su) — 0; 

dv du^ dv^ du dx dy dt 

a a a a 

-auv zr- (St) + au^ (Sv) + av^ (Sv) - av ^ (St) - (Al-69) 
dx du dv dt 

, a a^ a^ as a a a 

-av^ zr- (St) -aSv + kz^ (Sv) + A: ^ (Sv) -2k ^-u^^ (Sv) - v ^ (Sv) - ^ (^v) = 0; 
dy du^ dv^ dv dx dy dt 

a a a a 

-ku zr- (St) -kvz^ (St) + 2kz^ (Su) -kz^ (St) 0; (Al-70) 
dx dy du dt 

a a a a 

-ku zr- (St) -kvz^ (St) + 2kz^ (Sv) -kz^ (St) = 0; (Al-71) 
dx dy dv dt 

a a 

2A;^(^m)-i-2A:^(^v)=:0; (Al-72) 
dv du 

a dB d dB d 

-2au -zr- (St) -au^-2av z^ (St) -av ^-2az- (St) - (Al-73) 
dx du dy dv dt 

_^d^B _^d^B _^^dB _^^dB _^dB _^ 
du^ dv'^ dx dy dt 
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a 




3 




9 2 


3 3 


3 




8u — -{uv 




m- 


dx 


{Sx) + u 




^iSt)-v^iSx)- 
ax ay 


3? 


{Sx)y, 




d 




d 




3 


3 2 3 


3 




5v — —{uv 


dx 


{5t)- 


dx 


{Sy)-v 


dy 




dt 


{Sy}). 



3A 3A 3^A 3^A 3A 3A 3A 

-au -z: av ^ 2aA - k -z—r- - -k -z-^ + Mrr— + vrr— + — =0. (Al-74) 

au av au^ av^ dx dy dt 

Equation (Al-74) is simply Fokker - Planck equation for A. 
We solve (Al-66 - Al-67) and find 5u, 5v 

7i 7\ 7\ 7\ 

(Al-76) 
(Al-77) 

IX ay oi ay ai 

This gives for (Al-67 - Al-73) 

3 3 3 3 3 

-lauv ^ {5t) -au^ {St) - 2au^ ^ (St) -a^ {5x) -2k ^ (St) - (Al-78) 
dy dt dx dt dx 

dB 3^ 32 2 32 3^ 32 

-2k ^ + 2uv (St) - 2uv z^-^ (dx) + uv^ ^ (St) + u — (St) - 2u (Sx) + 
du dtdy dxdy dy^ dt^ dtdx 

g2 g2 g2 g2 g2 g2 ^2 

+2u'^v .-^ (St) + 2u^ (St) -u^^ (Sx) + u^^ (St) - 2v (Sx) - — (Sx) - — r (Sx) - 0; 
dxdy dtdx dx^ dx^ dtdy dy^ dt^ 

3 3 3 3 3 

-2auv ^ (St) -av — (St) - 2av^ ^ (St) -a^ (Sy) -2k ^ (St) - (Al-79) 
dx dt dy dt dy 

dB 32 32 , 32 32 

-2k rr- + 2uv (St) - 2uv z-z— (Sy) + 2uv^ z-z— (St) - 2u -z—- (Sy) + 
dv dtdx dxdy dxdy dtdx 

3 3 3 3 

-5te — (^0 -Zkv — (St) + 2k — (Sx) -3k — (St) ^ 0; (Al-80) 
dx dy dx dt 

3 3 3 3 

-3ku — (St) -5kv — (St) + 2k — (Sy) -3k — (St) 0; (Al-81) 
dx dy dy dt 

3 3 3 3 

-2ku — (St) - 2kv — (St) + 2k — (Sx) + 2k — (Sy) = 0; (Al-82) 
dy dx dy dx 

3 3B 3 3B 3 

-2au — (St) - au — 2av — (St) - av — 2a — (St) - (Al-83) 

dx du dy dv dt 

d^B d^B dB dB dB 

-k ^r-^ - k -z—r- + M — + V — + — =0. 

3m^ 3v^ dx dy dt 

Now we can collect similar terms in (Al-80 - Al-82) and so split them into nine equations: 

3 

-5k z- (St) = 0; (Al-84) 
dx 

d 

-3k — (St) = 0; (Al-85) 
dy 
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d d 

2k ^ (Sx) -3k — (St) = 0; (Al-86) 
dx at 

a 

-3k ^ {St) := 0; (Al-87) 
ax 

a 

-5A: ^ (^0 = 0; (Al-88) 
dy 

3 a 

2k — (Sy) -3k — (St) = 0; (Al-89) 
ay at 

a 

-2k ^ (St) = 0; (Al-90) 
oy 

a 

-2k ^ (St) = 0; (Al-91) 
ax 

a a 

2k^iSx) + 2k^iSy)^0. (Al-92) 
ay ax 

From (Al-84 - Al-85), (Al-87 - Al-88), (Al-90 - Al-91) we see, that St = St{t), which results in 
ftirther simplifications 

a a dB a^ a^ 

-au {St) - a {Sx) - 2yfc — - 2mv {Sx) + {St) - 2u :-— {Sx) - (Al-93) 
at at au axay at^ atax 

& ^'^^ - My ^'^^ - W ^'"^ " ^ ^'"^ = 

-av ^ {St) - a ^ {Sy) -2k^-2uv^ {Sy) - 2m ^ {Sy) - (Al-94) 
at at av axay atax 

^ ^'y^ ^ ' w ^''^ - &y ^'y^ - ^ ^'y^ - 1 ^'y^ = «^ 

a a 

2k ^ {Sx) -3k ^ {St) = 0; (Al-95) 
ax at 

a a 

2A: ^ (Jy) -3k ^ {St) = 0; (Al-96) 
ay at 

a a 

2)fc::-(Jx) + 2yfc::— (^y) = 0; (Al-97) 
ay ax 

dB dB a d^B d^B dB dB dB 

-au av- 2a — {St) - k ^—r - k ^—r + u — + v — + — =0. (Al-98) 

au av at au^ av^ ax ay at 



We integrate (Al-95 - Al-96) and find 



a 

Sx = C + 3/2x — {Sty, (Al-99) 
at 

a 

Sy^D + 3/2y ^ {St); (Al-lOO) 
at 
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where C = Ciy, t),D- D{x, t). We substitute these expressions to (Al-93 - Al-94), (Al-97 - Al-98) and 
obtain 

^ 3 3C dB 3^ 

-ZIlax ^ (St) -au^ (St) -a — -2k — - 3/2x ^ (St) - (Al-101) 
at'^ dt dt du 3?^ 

^ 3^ ^ 32c 232c d^C ^ 

-2u ^ {5t)-2v - V ITT - ITT = 0; 
3^2 3f3y 3y2 dt^ 

32 3 3D dB 33 

-3l2ay — (St) - av — (St) -a — -2k — - 3/2y — (St) - (Al-102) 
dt-^ dt dt dv dt' 

-2u - — - 2v ^ (St) - — = 0; 
3?3x 3x2 

3C dD 

2k^ + 2k^^0; (Al-103) 
dy dx 

dB dB 3 325 32^ dB dB dB 

-au av 2a ^ (St) - k ^-^ - k -i-m — + v — + — =0. (Al-104) 

du dv dt du^ dv^ dx dy dt 

We find ^ from (Al-101) and ^ from (Al-102): 
du dv 

35 1 32 3 dC 33 

— = - (-3/4ax ^ (St) - ll2au ^ (St) - \l2a ^ - 31 Ax ^ (St) - (Al-105) 
du k dt-^ dt dt dt"* 

32 32c ,32c 32c 



dB 1 , _ , . 32 _ , , 3 _ , . dD _ , . 3 
37 



,3 



- (-3/4ay ^ (St) - l/2av ^ (St) - ma ^ - 31 Ay —r (St) - (Al-106) 
k dt^ dt dt dt'' 



d^D jd^D 32 d^D 

-'^did-x-'"-' d^-^dT^ W'- 



Differentiating (Al-105) by v we have 



32^ _ 1 ^ d^C d^C ^ ^^^^ 
3m3v k 3y2 dtdy 

Differentiating (Al-106) by u we have 

325 1 d^D d^D 

^ = T^-''l^-^^- (Al-108) 
dudv k dx-^ dtdx 

We know, that C = C(y, t), D = D(x, t) and so we conclude from (Al-107 - Al-108) 

C = Ciy + E; (Al-109) 



D^C2X + F; (Al-110) 



where E = E(t), F = F(t), €^+€2^ 0. 
We find derivatives of B. 



dB I 32 3 dE 33 32 32£ 

— - — (3/2ax — (St) + au — (St) + a — - 3l2x — (St) - 2m — (St) - — ); (Al-11 1) 

du 2k dt-^ dt dt dt' dt-^ dt-^ 
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d^B d^B 

= 0. (Al-113) 



dudv dvdu 
Integration of (Al-111 - Al-112) gives 

1 , a dE a^ , a^ b^e 

B^G+ — {3/2uax zr^ (St) + lllau^ ^ {6t) + ua^- 3/2ux zr-r (St) - m ^ (St) - u ^r— ) + (Al-1 14) 
2k dt-^ at at at'' at^ at^ 

1 a^ , a dF a^ , a^ av 

^ 2^ a^ ^^^^ ^ 1^^^^ aF ^''^ ^ ar - ^^^^^ a7^ ^^^^ " ^ a^ ^^^^ - ^ a^^^ 

where G = G{x, y, t). 

Substitution of (Al-1 14) to (Al-98), collecting and equating to zero similar terms by u, v gives 

, , a^ , 1 a£ , a"* , a^s do 

3/4a^k-^x — (St) + l/2a^k-^ — - 3/4r'x — (St) - l/2r' — 4- — 0; (Al-115) 
at^ at at^ at^ ax 

9 r a , a^ 

l/2a^k~^ ^ (St) - 5/4r' ^ (St) ^ 0; (Al-116) 
at ot^ 

, , a^ , , aF , a'^ , av ao 

3/40^^;-^ ^ (^0 + l/2a^A;-' ^ - 3/4rV ^ (St) - l/2r' ^4-^=0; (Al-117) 
at^ at at^ at^ ay 

9 I a , a^ 

l/2a^A:-i — (St) - 5/4r' — (St) ^ 0; (Al-1 18) 

at at^ 

a a^ ac 

-5, (^0 + 2 ^(.0+3^^0; (Al-119) 

Integrate (Al-115), (Al-117) and obtain following expression (H - H(t)): 

G^H- - 3/8a^x2 — (St) + \l2xa^ - 3/8x^ ^ (^0 - l/2x — - (Al-120) 

1 f , , 32 ,dF ,34 av^ 

- I 3/8a^/ - (St) . mya- - - 3ISy^ - (St) - 1/2, — 



Substitution of (Al-120) to (Al-119), collecting and equating to zero terms by x, y gives 

d^E d^E 
-l/2a^k~^ — + l/2r' ^ = 0; (Al-121) 
at-^ at^ 

3^ 3^ 
-3/Sa^k-^ ^ (St) + 3/8r' ^ (St) 0; (Al-122) 
at' at-' 

, , 3V , 3V 

-l/2a^r^ ^ l/2r' ^ = 0; (Al-123) 
3^2 3?'^ 

3^ 3^ 
-3ma^k-^ — (St) + 3/8r' — (^0 = 0; (Al-124) 
at^ at^ 

3 3^ dH 

-9,(^0 + 2^(.0 + ^^0. (Al-125) 
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From (Al-116), (Al-118), (Al-122), (Al-124) we conclude, that 

St ^ const ^ Cj,. (Al-126) 

From (Al-121), (Al-123) we conclude, that 

E^C^ + C^t + C^e-''' + C^e'"■, (Al-127) 

F = Cg + Cgf + Cioe'"' + Cue"'; (Al-128) 
From (Al-125) and (Al-126) we see, that 

H = Ci2. (Al-129) 

We obtain using (Al-126 - Al-129) and backward substitution final expressions for variations: 

Sn = -l/lak'^unCs - l/lak'^nC^ - l/la^k'^xnCs - (Al-130) 

-l/2a^k~^ynC9 - a^k~\nC-je"' - a^k~\nCne"' + nCn + M 

5x ^ yCx + tC^ + Ci + C^e'"' + Cje"'; (Al-131) 

Sy = xC2 + tCg + Cg + Cioe""' + Cue"'; (Al-132) 

Su = -aCee'"' + aCjc"' + vCj -I- C5; (Al-133) 

Sv ^ -aCioe~"' + aCne"' + mCj + C9; (Al-134) 

St^C^. (Al-135) 

This ends calculations. 



